Structure and rotations of the Hoyle state 
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The Hoyle state is one of the most interesting and important challenges in nuclear physics. This 
excited state of the 12 C nucleus plays a key role in the fusion of three alpha particles to form carbon 
in red giant stars. In this letter we present ab initio lattice calculations which unravel the structure 
of this state and find evidence for a low-lying spin-2 rotational excitation. For the structure of the 
12 C ground state and first spin-2 state, we find a compact triangular configuration of alpha clusters. 
For the Hoyle state and second spin-2 state, we find a bent-arm or obtuse triangular configuration 
of alpha clusters. We calculate electromagnetic transition rates among the various low- lying states. 
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The carbon nucleus 12 C is produced by fusion of three 
alpha particles inside red giant stars. Without any res- 
onant enhancement, however, the triple- alpha reaction 
rate is too low to explain the abundance of carbon in our 
universe. In the early 1950's, Opik and Salpeter noted 
independently that the first step of merging two alpha 
particles is enhanced by the formation of 8 Be [T]-[3]. The 
ground state of 8 Be is a resonance with energy 92 keV 
above the 4 He- 4 He threshold with a width of 2.5 eV. A 
year later, however, Hoyle realized that this enhancement 
is still not enough to explain the observed carbon abun- 
dance. To resolve the discrepancy, Hoyle predicted a new 
positive-parity resonance of 12 C just above the combined 
masses of 8 Be and 4 He [4]. 

About three years later, Cook, Fowler, Lauritsen, and 
Lauritsen, experimentally observed a J n = + state 
278 keV above the 8 Be- 4 He threshold 0. This 0^ ex- 
cited state has a width of 8.5 eV, and is now commonly 
known as the Hoyle state. The triple- alpha reaction is 
completed when the Hoyle state decays electromagnet- 
ically to the 2+ state and then subsequently to the 0^ 
ground state. Around the time of the experimental con- 
firmation of the Hoyle state, Morinaga conjectured that 
the structure of excited alpha-nuclei such as the Hoyle 
state maybe be non-spherical [6] . This would imply low- 
lying even-parity rotational excitations. But there are 
also other ideas for the Hoyle state structure such as a 
diffuse trimer of alpha particles [7]. Recently the spin- 
2 Hoyle state has attracted considerable interest as the 
focus of several experimental investigations [SHE] • 

Recently we have presented an ab initio lattice calcu- 
lation of the Hoyle state [13] . In that work the low-lying 



spectrum of 12 C was calculated using the framework of 
chiral effective field theory and Monte Carlo lattice sim- 
ulations. In this letter we continue these investigations 
and present ab initio lattice results that answer some 
questions about the structure of the Hoyle state and the 
existence of rotational excitations. We find evidence for 
a low-lying spin-2 rotational excitation of the Hoyle state. 
For the Hoyle state and its spin-2 excitation, we find 
strong overlap with a bent-arm or obtuse triangular con- 
figuration of alpha clusters. This is in contrast with the 
12 C ground state and the first spin-2 state, where we see 
strong overlap with a compact triangular configuration 
of alpha clusters. We also calculate the electromagnetic 
transition rates among the low-lying even-parity states 
of 12 C. Our lattice results can be compared with other 
recent theoretical calculations for the low-lying spectrum 
of 12 C using the no-core shell model [I4j [15] and varia- 
tional calculations using Fermionic Molecular Dynamics 

maun- 

Chiral effective field theory treats the interactions of 
protons and neutrons as a systematic expansion in pow- 
ers of nucleon momenta and the pion mass. A recent 
review can be found in Ref. [18 . The low-energy ex- 
pansion is organized in powers of Q, where Q denotes 
the typical momentum of particles and is treated on the 
same footing as the mass of the pion. The most impor- 
tant contributions to the nuclear Hamiltonian come at 
leading order (LO) or O(Q ), while the next-to-leading 
order (NLO) terms are 0(Q 2 ). In the lattice calcula- 
tions presented here, we include all possible interactions 
up to next-to- next-to-leading order (NNLO) or 0(Q 3 ). 

In this analysis we use a periodic cubic lattice with lat- 
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tice spacing a = 1.97 fm and total length L = 12 fm. In 
the time direction, we use lattice time step a t = 1.32 fm 
and vary the propagation time L t to extrapolate to the 
limit L t — >• oo. The nucleons are treated as point-like 
particles on lattice sites, and interactions due to the ex- 
change of pions and multi-nucleon operators are gener- 
ated using auxiliary fields. Lattice effective field theory 
was originally used to calculate the many-body proper- 
ties of homogeneous nuclear and neutron matter p~9j [20] . 
Since then the properties of several atomic nuclei have 
been investigated [2TJ [22] . A recent review of the liter- 
ature can be found in Ref. [23] . 

Euclidean time propagation is used to project on to 
low-energy states of our interacting system. Let H be 
the Hamiltonian. For any initial quantum state the 
projection amplitude is defined as the expectation value 
(e~ Ht )y For large Euclidean time t, the exponential 
operator e~ Ht enhances the signal of low-energy states. 
Energies can be determined from the exponential decay 
of these projection amplitudes. The first few time steps 
and last few time steps are evaluated using a simpler 
Hamiltonian #su(4) based upon Wigner's SU(4) symme- 
try for protons and neutrons [24] . This Hamiltonian is 
computationally inexpensive and is used as a low-energy 
filter before starting the main calculation. This tech- 
nique is described in Ref. [23] . 

In Table [I] we present lattice results for the ground 
state energies of 4 He and 8 Be up to NNLO. The method 
of calculation is nearly the same as that described in 
Ref. [13, 22, 25 . The higher-order corrections are com- 
puted using perturbation theory. The coefficients of 
the nucleon-nucleon interactions are set by fitting to low- 
energy scattering data. In our calculations the NNLO 
corrections correspond with three-nucleon forces. A de- 
tailed description of the interactions at each order can 
be found in Ref. [25]. We have used the triton binding 
energy and the weak axial vector current to fix the low- 
energy constants cd and ce entering the three-nucleon 
interaction. 

In comparison with the calculations in Ref. [13 , some 
improvements have been made using higher-derivative 
lattice operators which eliminate the overbinding of the 
leading order action when calculating larger nuclei such 
as 16 0. The details of this improved action will be dis- 
cussed in a forthcoming publication. The error bars in 
Table I are one standard deviation estimates which in- 
clude both Monte Carlo statistical errors and uncertain- 
ties due to extrapolation at large Euclidean time. We 
see that the binding energy results for 4 He and 8 Be at 
NNLO are in agreement with experimental values. 

In our projection Monte Carlo calculations we use a 
larger class of initial and final states than considered in 
previous work. For the calculation of 4 He we use an 
initial state with four nucleons, each at zero momentum. 
For the calculation of 8 Be we use the same initial state 
as 4 He, but then apply creation operators after the first 



TABLE I: Lattice results and experimental values for the 
ground state energies of 4 He and 8 Be. All energies are in 
units of MeV. 
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FIG. 1: Lattice results for the 12 C spectrum at leading order. 
In Panel I we show results from three different initial states, A, 
B, and A, each approaching the ground state energy. In Panel 
II we show results starting from three other initial states, C, 
D, and A. These trace out an intermediate plateau at energy 
about 7 MeV above the ground state. 



time step to inject four more nucleons at zero momentum. 
The analogous process is done to extract four nucleons 
before the last step. This injection and extraction pro- 
cess of nucleons at zero momentum helps to eliminate 
directional biases caused by initial and final state mo- 
menta. 

We make use of many different initial and final states 
to probe the structure of the 12 C states. In all of the 12 C 
states investigated here we measure four-nucleon correla- 
tions by calculating the expectation value of p 4 , where p 
is the total nucleon density. We find strong four-nucleon 
correlations consistent with the formation of alpha clus- 
ters. In Fig. [I] we present lattice results for the energy 
of 12 C at leading order versus Euclidean projection time 
t. For each of the initial states A, B, C, and D, we 
start with delocalized nucleon standing waves and use a 
strong attractive interaction in .Hsu (4) to allow the nu- 
cleons to self-organize into a nucleus. For initial states 
A and A, we use alpha cluster wavefunctions to recover 
the same states found using initial states A, B, C, and D. 
For these calculations, the interaction in #su(4) is n °t as 
strong and the projected states retain their original alpha 
cluster character. 

In Panel I, we show results from three different initial 
states, A, B, and A, each approaching the ground state 
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FIG. 2: This shows initial state A, a wavefunction consisting 
of three alpha clusters formed by Gaussian packets centered 
on the vertices of a compact triangle. There are a total of 12 
equivalent orientations of this configuration. 




FIG. 3: This shows initial state A, a wavefunction consisting 
of three alpha clusters formed by Gaussian packets centered 
on the vertices of a bent-arm or obtuse triangular configura- 
tion. There are a total of 24 equivalent orientations of this 
configuration. 



energy, —96(2) MeV. For initial state A, we start with 
four nucleons each at zero momentum, apply creation 
operators after the first time step to inject four more 
nucleons at rest, and then inject four more nucleons at 
rest after the second time step. The reverse process 
is used to extract nucleons for final state A. The same 
scheme is used for initial state B, though the interactions 
in .Hsu (4) used are not as strongly attractive as those for 
A. 

For initial state A, we use a wavefunction consisting 
of three alpha clusters as shown in Fig. [2] The alpha 
clusters are formed by Gaussian packets centered on the 
vertices of a compact triangle. In order to construct 
eigenstates of total momentum and lattice cubic rota- 
tions, we consider all possible translations and rotations 
of the initial state. There are a total of 12 equivalent 
orientations of this configuration. We do not find fast 
convergence to the ground state when starting from any 
other configuration of alpha clusters. From this we con- 
clude that the alpha cluster configurations in Fig. [2] have 
the strongest overlap with the Of ground state of 12 C. 
The fact that it is an isosceles right triangle rather than 
an equilateral triangle is just an artifact of the lattice 
spacing. 

In Panel II of Fig. [I] we show leading-order energies 
for three different initial states, C, D, and A, each ap- 
proaching an intermediate plateau at —89(2) MeV. If 
Euclidean time is taken to infinity, these curves eventu- 
ally approach the ground state energy like the curves in 
Panel I. However it is clear that a different state is first 
being formed which is not the ground state. We identify 
the + state in this plateau region as the 0^ Hoyle state. 
The common thread connecting each of the initial states 
C, D, and A, is that each produces a state which has an 
extended or prolate geometry. This is in contrast to the 
oblate triangular configuration in Fig. [2] 

For initial state C, we take four nucleons at rest, four 
with momenta (27r/L, 27r/L, 27r/L), and four with mo- 
menta (— 27r/L, — 27r/L, — 2tt/L). For initial state D, 
we use a similar configuration with four at rest, four 
with momenta (2ir/L, 2tt/L, 0), and four with momenta 
(— 27r/L, — 27r/L, 0). For initial state A, we use a set 
three alpha clusters formed by Gaussian packets centered 
on the vertices of a bent-arm or obtuse triangular con- 



TABLE II: Lattice results for the low- lying even-parity states 
of 12 C compared with the experimental results in units of 
MeV. 





of 


2+(£+) 


o 2 + 


2f(£ + ) 


LO [O(Q )] 


-96(2) 


-94(2) 


-89(2) 


-88(2) 


NLO [0{Q 2 )\ 


-77(3) 


-74(3) 


-72(3) 


-70(3) 


NNLO [0(Q 3 )] 


-92(3) 


-89(3) 


-85(3) 


-83(3) 


Experiment 


-92.16 


-87.72 


-84.51 


-82.6(1) [EJCEo] 
-82.32(6) [IT] 
-81.1(3) [9] 



figuration as shown in Fig. [3J There are a total of 24 
equivalent orientations of this configuration. We do not 
find the same plateau starting from other configurations 
of alpha clusters. We conclude that the configurations 
in Fig. have the strongest overlap with the Of Hoyle 
state of^ 2 C. 

We use the same multi-channel method developed in 
Ref. [13 to find a spin-2 excitation above the ground 
state as well as a spin-2 excitation above the Hoyle state. 
In both cases we are taking the E + representation of the 
cubic rotation group on the lattice. We show the re- 
sults for the binding energies of the low-lying even-parity 
states of 12 C in Table [n] We find that the binding ener- 
gies at NNLO are in agreement with experimental values. 



In Table III we present results at leading order for the 
root-mean-square charge radius and quadrupole moment 
of the even-parity states of 12 C. We also show experi- 
mental values where available. In this study we compute 
electromagnetic moments only at leading order. We note 
that moments such as the charge radius for resonances 
above threshold are dependent on boundary conditions 
used to regulate the continuum-state asymptotics of the 
wavefunction. We avoid this problem because all of the 
low-lying states are bound at leading order. One expects 
that as the higher-order corrections push the binding en- 
ergies closer to the triple alpha threshold, the correspond- 
ing radii will increase accordingly. A detailed study of 
these resonances as a function of finite volume size will 
be investigated in future work. We find good agreement 
with the experimental value for the 2\ quadrupole mo- 
ment. The difference in signs for the electric quadrupole 
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TABLE III: Lattice results at leading order and experimen- 
tal values for the root-mean-square charge radius and the 
quadrupole moments for 12 C. 





LO 


Experiment 


r(0+) [fm] 


2.2(2) 


2.47(2) [26] 


r(2+) [fm] 


2.2(2) 




Q(2+) [efm 2 ] 


6(2) 


6(3) [27] 


r(0+) [fm] 


2.4(2) 




r(2+) [fm] 


2.4(2) 




Q(2+) [efm 2 ] 


-7(2) 





TABLE IV: Lattice results at leading order and experimen- 
tal values for electromagnetic transitions involving the even- 
parity states of 12 C. 





LO 


Experiment 


B(E2,2+ ->0+) [e 2 fm 4 ] 


5(2) 


7.6(4) [29] 


B{E2,2+ -S-0+) [e 2 fm 4 ] 


1.5(7) 


2.6(4) [29] 


B(E2,2+ ->0+) [e 2 fm 4 ] 


2(1) 




B(£2,2+^0+) [e 2 fm 4 ] 


6(2) 




m(£0,0j ->• 4 ") [e fm 2 ] 


3(1) 


5.5(1) na 



moments of the two spin-2 states reflects the oblate shape 
of the 2+ state and prolate shape of the 2^ state. 

The leading order results for the electromagnetic tran- 
sitions among the even-parity states of 12 C are shown 
in Table IV The definitions for these quantities can be 
found in Ref. [28]. The agreement with available ex- 
perimental values is reasonable. The lattice results at 
leading order have a tendency to be somewhat smaller 
than experimental values. This presumably reflects the 
greater binding energies and smaller radii of the nuclei at 
leading order. We also predict electromagnetic decays 
involving the 2^ state that may be measured experimen- 
tally in the near future. 

In summary we have presented ab initio lattice cal- 
culations which show the structure of the Hoyle state 
and find evidence for a low-lying spin-2 rotational ex- 
citation. For the ground state and first spin-2 state, 
we find mostly a compact triangular configuration of al- 
pha clusters. For the Hoyle state and second spin-2 
state, we find a bent-arm or obtuse triangular config- 
uration of alpha clusters. We have calculated charge 
radii, quadrupole moments, and electromagnetic transi- 
tions among the low-lying even-parity states of 12 C at 
leading order. All of the results are in reasonable agree- 
ment with experimental values. More work is still needed 
such as calculations using smaller lattice spacings. How- 
ever these results provide a deeper understanding of the 
structure and rotations of the Hoyle state starting from 
first principles. 
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